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Abstract 

The 0(3) nonlinear a model is studied in the disordered phase, using the 

techniques of the effective action and finite temperature field theory. The 

nonlinear constraint is implemented through a Lagrange multiplier. The finite 

temperature effective potential for this multiplier is calculated at one loop. 

The existence of a nontrivial minimum for this potential is the signal of a 

disordered phase in which the lowest excited state is a massive triplet. The 

mass gap is easily calculated as a function of temperature in dimensions 1, 

2 and 3. In dimension 1, this gap is known as the Haldane gap, and its 

temperature dependence is compared with experimental results. 
75.10.Jm, ll.10.Lm 
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The 0(3) nonlinear a model (NLcr) offers a good description of quantum antiferromag- 
nets at long wavelengths, as suggested by renormalization group analysesS. The mapping 
between the nonlinear a model and the Heisenberg model was first established by Haldanei, 
who also argued that a topological Hopf term should be present in the case of half-integer 
spins in one dimension! (this distinction between half-integer and integer spins disappears 
in two dimensions). For integer spins, in which case the Hopf term is absent, Haldane ar- 
gued that the lowest excitations should exhibit a mass gap. Accordingly, these excitations 
(the magnons) are no longer Goldstone modes and the rotation symmetry of the spins is no 
longer spontaneously broken: the symmetry has been dynamically restored by quantum fluc- 
tuations. This phenomenon has found an experimental realization in quasi one-dimensional 
antiferromagnets in their disordered phasei. 

In this paper we calculate the temperature dependence of this mass gap in dimensions 
1 to 3 within the framework of the finite temperature effective action of the nonlinear a 
model. The use of this technique is facilitated by the Lorentz invariance of the NLct model, 
which ensures that the magnon speed is not affected by quantum corrections. 

The 0(3) nonlinear a model describes the dynamics of a unit vector (p representing the 
local direction of the staggered magnetization. The Euclidian action obtained from the 
Heisenberg model is 

S = Yg I dTddx {^ {dT0)2 + {dl0)2 } (1) 

wherein v is the spin-wave velocity and g a coupling constant, the constraint ip 2 = 1 
is understood. If we choose a set of units in which the speed v is unity, implement the 
constraint (p 2 = 1 by a Lagrange multiplier a, and rescale the field p by y/g in order to bring 
its kinetic term into the standard form, we obtain the equivalent action 

S = J drd d x {\d^pd^-\a{p 2 - l/g)} (2) 

wherein the index /i goes from to d. The spin-wave velocity v may always be restored later 
using dimensional analysis. 



The action (|2]) is symmetric with respect to global rotations of the field (p, which is 
the relevant order parameter. The degenerate classical ground states are configurations 
in which the field <p is uniform. Such ground states break the full symmetry group 0(3) 
down to 0(2). Goldstone's theorem would then imply the existence of two massless modes, 
corresponding to the two magnon polarizations. However, the quantum ground state, to 
which Goldstone's theorem refers, does not break the 0(3) symmetry if we are in a disordered 
phase. The symmetry of the quantum ground state may be determined by calculating the 
one-loop effective potential V(a) for the field a, following the discussion by Coleman!. If 
this potential has a global minimum at some nonzero value <To of the field a, then we must 
expand the field a around its ground state value cr . The latter then gives rise to a mass 
term ^cr o 2 in the action (|2|). As a consequence, the lowest lying excitations are no longer a 
pair of massless Goldstone modes (the infrared instability has been removed) but a massive 
triplet: the elementary excitations of the three-component field (p with mass gap A = y/aQ. 

The one-loop effective potential V(a) is obtained by summing an infinite sequence of 
Feynman graphs, each with n insertions of a and n internal (p lines. This amounts to evalu- 
ating the partition function with a constant, which is done simply by Gaussian integration. 
The result of this standard calculation is 



where we have included the tree- level potential cr/2g (we use the convention h = 1 through- 
out). In order to find out whether V(a) has a minimum, we consider instead the derivative 



The extension of this formula to a finite temperature T is obtained by replacing the inte- 
gration over pq by a sum over Matsubara frequencies u n = 2imT: 




(3) 






(5) 



Using the standard result ((3 = 1/T) 



3 



n ™n 



we may write 



1 



At zero temperature, the above reduces to 



2V>) = l_ 3Kd f J^L (8) 
g Jo VP + °" 



where K^ 1 = 2vr and Kj 1 = 2 d ~ V /2 r(d/2) if d > 1. An ultraviolet cutoff A has been 
introduced. We point out that the NLa model is a renormalizable field theory only in d = 1; 
this means that our treatment for d > 1 will only be valid in the limit of large distances. 
Without taking A to infinity, we will consistently neglect terms of order A -1 or more. This 
being understood, finite temperatures will not introduce any new dependence on A: Indeed, 
the difference Vs = V(a, T) — V(a, 0) is given by 

2 y» = -m d J^^M= {coth (Vy^M-T) - l} (9) 

Since the function cothr — 1 decreases exponentially^ as x — > oo, the above integral is UV 
finite and the cut-off may be removed. The dependence on A is contained solely in V(a, 0). 
A change of integration variables yields 



2V» = -T d - l Uy) (10) 



where y = ■sfo /T and 



a-its r A „ coth(|^) - 1 



We find that 



Uv) = W- 1 K d l dz (z2 X _ 2 l)i d/ 2 (11) 



q OO 

f 1 (y) = -Y, K o(ny) (12a) 

f 2 (y) = -{b~ log(2sinh(y/2))} (12b) 

/3(») = jE-*) ( 12c ) 

" n=l " 
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(Ko and K\ are modified Bessel functions.) 

Let us now look in more detail at specific cases, starting with d — 1. This has been 
done explicitly by Colemani for the CP 1 model, which is equivalent to the 0(3) nonlinear 
a model. From Eq. (|8|) we obtain, at zero temperature,0 

2V'(a) = - + A log (a/A 2 ) (d = 1) (13) 

We see that the equation V'(cr) = always has a solution, at the point ao = A 2 exp — (2ir/3g), 
which implies a zero temperature Haldane gap Ao = A exp —ir/3g. 
From Eqs. fllCf ) and (0), the finite temperature result is then 

1 Q o OO 

2V'(a) = - + f- log (a/A 2 ) - - E K (n^/T) (14) 

9 27r V ' ^=1 

The finite temperature Haldane gap A = ^/cfo is therefore obtained through the solution of 
the equation 

oo 

ilog(5) = E^oMA) (15) 

n=l 

where we have defined the reduced gap 5 = A/A and the reduced temperature t = T/A . 
This gap equation may be solved numerically. The solution is shown on Fig.l. The gap is 
seen to maintain its T = value until about T ~ 0.2A , after which it takes off linearly. 
We naturally expect the gap to grow with temperature since A, being the inverse of the 
correlation length, is in some sense a measure of disorder. This temperature dependence 
seems confirmed by experiments!. 

The dependence of A on T is not qualitatively affected by the introduction of anisotropy 
in the a model. Easy-axis anisotropy is introduced via a mass term |M 2 <£>^ in the action 
(0), while easy-plane anisotropy is rendered by a mass term |M 2 (<^ 2 + <^ 2 ). The above 
calculation of the Haldane gap is easily adapted to the anisotropic case by shifting a by 
a for one (two) out of three components of (p in the easy-axis (easy-plane) case. The gap 
equation (|T5| ) is thus modified into 

y + 5 2 



\ log I ( J * 4 | = E K *H 52 + » 2 /t) + 2 E KoH/t) (16) 



n=l n=l 



where we have defined the reduced anisotropy \i = M/A and considered the easy-axis 
case. Even if the value of A is affected by the anisotropy, everything else being fixed, we 
have checked that the t dependence of the reduced gap 5 is not changed appreciably by the 
introduction of an anisotropy even as large as fi = 1. However, an easy-axis anisotropy splits 
the Haldane triplet into a doublet with gap A x = A and a sing let with gap A 2 = VA 2 + M 2 . 
The weighted average of the two gaps, normalized to its zero temperature value, is shown 
on the lower curve of Fig.l for \i = 1.9. This value of /i, appropriate for NENP, is obtained 
from the observed values of Ai (14K) and A 2 (3OK)0. 

We next consider the two-dimensional case. At T = this was studied by Deser and 
Redlichl for the CP 1 model. The result of the integration (|) is 

2 ^>4-i + ^ (rf = 2) (17) 

The interesting feature here is the appearance of a critical coupling g c = 27r/A below which 
the gap equation V'(cr) = has no solution. This is the signal of a phase transition in cou- 
pling space: below g c , the 0(3) symmetry is spontaneously broken and massless Goldstone 
modes appear. At finite temperatures we have, following Eq. (0): 

1 1 3T r i 

2V'(a) = + — log 2sinh(y5 r /2T) (18) 

9 9c rc l 

The gap equation has now a solution for all nonzero values of T, indicating that the ordered 
phase does not survive at finite temperatures, as is already known from RG calculations^. 
Let us define the difference r] = l/g c — 1/ g. The solution of the gap equation is then 

A = 2Tsinh~ 1 [| exp(vrr//3T)] (19) 

If r] > (disordered phase) the gap tends towards the constant 27rrj/3 as T — >• 0, i.e., r\ = 
3A /27r. If, on the contrary, rj < (ordered phase) the gap behaves like A w Texp — ir\r]\/3T 
as T — > 0, i.e., it vanishes faster than any power law. At the critical value rj = the gap 
vanishes linearly as T is lowered: 

A = 2Tsinh~ 1 (i) w 0.962 T (20) 
6 



These results are illustrated in Fig. 2. Similar results were obtained within the RG analysis 
of Ref.[ |1| in terms of the correlation length £ = 1/A. Experimental results are available in 



knowledge, the disordered phase (77 > 0) has not been observed yet. Although it is accepted 
that this phase does not exist at T = for the Heisenberg model with spin s > 1, nothing 



Let us point out that the cross-over region, i.e., the temperature range over which thermal 
fluctuations overtake quantum fluctuations, is characterized by the condition A ~ T which, 
on the phase diagram (g,T), is met on the curve T ~ 25. 3rj, according to Eq. (0). When 
T goes beyond this value the states over the gap are thermally excited or, stated differently, 
the correlation length £ exceeds the extent j3 of the Euclidian time direction; the fluctuations 
then happen mostly in the spatial directions. 

Let us finally discuss the three-dimensional case, first at T = 0. The integration @ 
yields 



Again, we have dropped terms of order 1/A. It is understood that a is much smaller than 
the cutoff A, so that V'(a) is increasing monotonically with o. Like in dimension 2, there 
is a critical coupling g c = 4-7r 2 /3A 2 below which the symmetry is spontaneously broken at 
T = 0. The difference with d = 2 is that this ordered phase survives at finite temperature. 
Indeed, following Eqs. ([l|) and fl2T|), the gap equation is 



The function f$(y) is monotonically decreasing, rather like an exponential, starting with 
/ 3 (0) ~ 0.5. If the above equation has no solution at T — (g < g c ), then an increase 
in temperature will bring a solution at some point, since the r.h.s. can only decrease as T 
increases. The transition to the ordered phase occurs when A = 0; the relation between g 
and T that follows is 



the ordered phase (rj < 0) and compare well with the theory (see Ref. [ |IJ). To the author's 



certain can be said about the spin-^ case. 






(22) 
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--- = T 2 / 3 (0) (23) 
9 9c 

We see that no matter how small g may be below g c , a critical temperature exists beyond 
which we recover a disordered phase. However, this result is not to be trusted at arbitrary 
low g, since we tacitly assumed that T<A. 

If we place ourselves above g c , denote by A the zero temperature gap and define reduced 
variables 5 = A/A , t = T/A and A = log(A/A ) + |(log2 — 1), the gap equation becomes 

Att 

A{ 1 _ _ f i og5 + —t 2 f 3 (5/t) = (24) 

Its solution for 5 as a function of t may be found numerically, and is shown on Fig. 3 
for two values of A. Note that we cannot rid ourselves of this extra parameter, but its 
logarithmic dependence on the cut-off causes only a mild non universality. At the critical 
value g = g c , we find that the gap vanishes linearly with T, like in dimension 2, albeit with 
mild corrections. 

In conclusion, calculating the one-loop effective potential for the Lagrange multiplier 
a yields the Haldane gap as a function of temperature in a rather direct way. The solu- 
tion is compatible with experimental results. The extension to two and three dimensions is 
straightforward although, to the author's knowledge, experimental results in the quantum 
disordered regime are then lacking. Since no rigorous argument has yet been found that pre- 
vents the existence of a quantum disordered phase in the two-dimensional spin-| Heisenberg 
model, the question remains open. 
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FIGURES 

FIG. 1. Magnitude of the reduced Haldane gap 5 = A/Ao as a function of t = T/Ao in d = 1 
(upper curve). Experimental results from Ref.[6] are also shown for CsNiCi3 and NENP (error 
bars are not included). The lower curve is the weighted average of the two gaps in the presence of 
a fi = 1.9 easy-axis anisotropy, appropriate for NENP. 

FIG. 2. Magnitude of the mass gap A as a function of temperature for d = 2 and three values 
(indicated) of rj. 

FIG. 3. Magnitude of the reduced mass gap 5 = A/Ao as a function of t = T/Aq in d = 3 for 
A = 4 (higher curve) and A = 5 (lower curve). 
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